Dissipation and Tunnelling in Quantum Hall Bilayers 
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We discuss the interplay between transport and intrinsic dissipation in quantum Hall bilayers, 
within the framework of a simple thought experiment. We compute, for the first time, quantum 
corrections to the semiclassical dynamics of this system. This allows us to re-interpret tunnelling 
measurements on these systems. We find a strong peak in the zero-temperature tunnelling current 
that arises from the decay of Josephson-like oscillations into incoherent charge fluctuations. In the 
presence of an in-plane field, resonances in the tunnelling current develop an asymmetric lineshape. 
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Transport in quantum Hall bilayers has been the sub- 
ject of much recent interest H&ElSS EH3 The bilay- 
ers consist of two closely-spaced parallel two-dimensional 
electron layers in a double quantum well. If the Lan- 
dau level fillings of the layers are v\ = = 1/2, then 
Coulomb interactions between the layers drive a transi- 
tion to a ground state in which the bilayer as a whole 
exhibits the quantum Hall effect [lj. This ground state is 
believed to have a broken U(l) symmetry 0; there is a 
macroscopically coherent phase associated with the elec- 
trons' layer degree of freedom. The ground state can be 
viewed as an easy-plane ferromagnet |3j or, equivalently, 
as an excitonic condensate^, E| . There are also analogies 
with Josephson junctionsyj. 

A series of remarkable experiments [3. Fl llfij have been 
used to probe the internal degrees of freedom of this 
strongly correlated system. They show evidence for in- 
terlayer coherence and the linearly dispersing Goldstone 
mode resulting from the broken U(l) symmetry 0, and 
also, most recently, "excitonic" superfluidity Q. 

There still remain questions concerning the interlayer 
tunnelling spectrum of the bilayers. One prominent fea- 
ture in the IV characteristic is a sharp peak in the 
tunnelling current for small biases (between 10 and 100 
/iV) 2]. At low temperatures, reducing the bias leads 
to a sharp rise in the tunnelling current. This increase 
is cut off below 10 /iV so that the current falls to zero 
at zero bias. Existing theories for interlayer tunnelling 
are restricted to the classical limit of the underlying spin 
model, and do not produce this feature at the low tem- 
peratures relevant to experiment (ksT ~ 2/ieV). 

In this paper we study interlayer tunnelling within the 
framework of a simple "thought experiment" . We fol- 
low the relaxation of an initial charge imbalance across 
the bilayer. Including quantum (1/S) corrections to the 
dynamics of the pseudospin model of the bilayer 0, we 
see that electron tunnelling across the bilayer generates 
density waves. This quantum dissipative process leads 
to a zero-temperature tunnelling current of the form 
I oc 1/V for a bias voltage V above a threshold Vb; this 
is consistent with experimental measurements. The long- 



wavelength density fluctuations in the bilayer have an en- 
ergy gap A sw ; in the absence of disorder the threshold 
Vq corresponds to an energy of this order. Introducing 
an in-plane magnetic field to our calculations causes the 
small bias feature to develop into an asymmetric resonant 
peak in the tunnelling current, similar to that reported 
in 0. Below the threshold Vq, we find that the intrinsic 
dissipation causes any macroscopically coherent charge 
oscillations to decay in time, consistent with the absence 
of Josephson-like oscillations in experiments. 

We work in the pseudospin picture of the bilayer 3]. 
The charge imbalance on the bilayer is given by the z- 
component of the magnetisation of the pseudospins. The 
system is a ferromagnet due to Coulomb exchange. Since 
a z component in the magnetisation incurs a capacitative 
energy cost, the ferromagnet has easy-plane anisotropy. 

We use a model of discrete spins on a lattice, begin- 
ning with a spin- 1/2 system in which each spin repre- 
sents a local single-particle state within the lowest Lan- 
dau level. We work with a large-5* version of this system. 
This large-S* generalisation may be treated as a coarse- 
graining procedure. The limit of large S is the classical 
limit for the spin system. We work with the Hamiltonian: 
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where Si — Sfhi = S(mf,rnf,rnf) is the spin operator 
on site i of a square lattice with spacing Cq = ^/2ttIb 
where Is = (hc/eB) 1 / 2 is the magnetic length. The ex- 
change pe and the strength of the on-site repulsion D 
were derived from microscopic considerations by Moon 
et aL0. The tunnelling between the layers enters the 
problem through A§as: the splitting of the "bonding" 
and "anti-bonding" single-particle states in the double 
well. In our thought experiment, a gate is used to con- 
trol the charge imbalance on the bilayer: this adds a term 
Hy = —SVj2i m i to the Hamiltonian. Typical values 
for the model parameters in physical bilayer systems are 
l B ^ 20nm, A S as ^ 90^K, p E ^ 0.5K, D ~ 30K. 

We consider the properties of the ferromagnetic phase 
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FIG. 1: Mean field (classical) spin trajectories are equipoten- 
tials on the surface of the spin sphere. 



of this model, which is believed to describe the experi- 
mental systems. In particular, the observation of a lin- 
early dispersing peak in the tunnelling conductance 
indicates that the experimentally accessible state has fer- 
romagnetic order — the peak arises from the Goldstonc 
mode of the system. We therefore ignore the quantum 
disordered phase that exists for D/peS 2 3> 1 

We will now outline the dynamics of the classical fer- 
romagnetic system before we discuss the quantum effects 
which are responsible for the existence of a dissipative 
tunnelling current. In the classical system, a spatially 
uniform spin configuration will remain uniform forever. 
The total magnetisation precesses along a trajectory of 
constant energy. We sketch the trajectories in Fig. for 
Asas D which is the physically relevant regime. 

Near the (x-polarised) ground state, the precession fre- 
quency is given by A sw /?i where A sw = [A S as(Asas + 
D)] 1 / 2 is the energy gap for spinwave excitations (den- 
sity waves of charge imbalance across the bilayer). As 
the magnetisation precesses along these trajectories, the 
charge imbalance on the bilayer oscillates around zero. 
Far away from the ground state, the spin precesses 
around one of the two maximal energy states which lie 
close to the S z axis. This yields a Josephson-like alternat- 
ing current / ~ eAsAS cos{eVt/Ti) where V is the voltage 
across the bilayer due to capacitative charging. We stress 
that this is valid only for large charge imbalance (large 
V). 

Trajectories through the saddle point at m = (— 1, 0, 0) 
mark the boundary between oscillations around the 
ground state and oscillations around the maximal energy 
states. The saddle point trajectory crosses the zz-plane 
at three points, as shown in Fig. ^ The angle 9q satis- 
fies cos 6*o = 1 — 2(Asas/-D). It corresponds to a voltage 
difference of Vb = 2A sw /e across the layers. 

In our thought experiment, we imagine using a gate to 
induce a uniform charge imbalance on the bilayer. In the 
spin picture, the magnetisation is tilted out of the easy 
(xy) plane. The bias is then instantaneously removed 
and the bilayer finds itself in a highly excited state. 

The resulting behaviour depends crucially on whether 
the initial gate voltage V is above or below the saddle 
point value Vq. We treat the two regimes in separate 



perturbation theories. In both calculations, we expand 
around the classical limit using spinwave theory in a 1/S 
expansion. The leading terms in this expansion result in 
the leading terms in the I-V relation. 

We begin by considering the situation for initial ener- 
gies above the saddle point (V > Vq). When the gate 
voltage is removed, the charge imbalance results in a po- 
tential difference across the bilayer equal to the original 
gate voltage V. As discussed above, this causes an alter- 
nating current with a frequency (eV/Ti). In this regime, 
we can treat the tunnelling Asas perturbatively. 

The quantum spin system is distinct from the classical 
one in that spatially uniform oscillations do not persist — 
the tunnelling term breaks the global spin rotation sym- 
metry so that long-wavelength modes are no longer pro- 
tected from decay by Goldstone's theorem. The uniform 
mode decays by transferring energy into spin waves with 
finite wavevectors. The magnetisation falls to a lower 
trajectory on the spin sphere, corresponding to a net 
transfer of charge across the bilayer. We can therefore 
obtain the d.c. tunnelling current by calculating the rate 
of this dissipative process. 

It is straightforward to derive a bosonic spinwave the- 
ory for our model. We use the Holstein-Primakov repre- 
sentation: Sf + iS" = (2S - a\aj) l l 2 aj, S z = S 



a]aj. 



Expanding around m z = (eV/D) = s'm9, the quadratic 
part of the Hamiltonian is then easily diagonalised in the 
Fourier basis to give 
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where the a q are bosonic operators describing the spin- 
wave modes. The spinwave dispersion is given by e q = 
[PE-f(q)(D + PE-f(q))} 1/2 where 7(g) = 4 - 2 cos^cn) - 
2 cos(<7j,cn). The dispersion is linear and gapless at small 
q. (Tunnelling should produce a small energy gap but 
this does not affect the perturbative calculation we de- 
scribe here.) The spinwave velocity is v — Ib^ttDpe) 1 / 2 ■ 
Observe that the q = mode has been singled out in the 
Hamiltonian, and its energy is not given by the long wave- 
length limit of £q. The quanta of this mode carry S z — 1 
while spin waves with finite wavevector have S z — 0. 

Starting from an initial state \i), the energy dissipation 
rate is given by 
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The dissipation arises from the destruction of one quan- 
tum in the q = mode, and the generation of multi- 
ple spinwaves during tunnelling across the bilayer. To 
leading order in l/S, a pair of spin waves with opposite 
momenta is excited (Fig. |2Jt) . The relevant vertex is: 
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FIG. 2: Decay processes for V > Vq. (a) dominant process for 
Q=0: Eq+E- q — eV . (b) Decay process at finite Q: eq = eV . 
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FIG. 3: Plots of the current I (arbitrary units) at Q — 
and at finite Q = 0.5{4irp B S/v) so that vQ ~ 150/iV. The 
two spinwave contribution (2sw) to the current gives rise to 
the asymmetric peak, as shown. The contributions at finite 
Q have been broadened by a Gaussian distribution for Q of 
width <jq = 0.1(47rp B .Syv). 



in which Asas = Asas exp(— sj Dj PeS 2 ) is a renor- 
malised tunnelling amplitude, and 72,^ = cos0[(u q + 
v q ) 2 + 2 sin 9 sec 2 6* — — v a ) 2 sec 4 6*1 where u a and zj 



are coherence factors: (u q + w^) 2 
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with u 2 — v 2 = 1. The ^-dependence arises because the 
x-component of the pseudospin depends on its angle with 
the easy plane, 9, as well as on the azimuthal angle <f>. 
This dependence is weak for 9 <C 1 when the charge im- 
balance is small compared to the Landau level filling. 

We now calculate the power dissipation T for a given 
initial voltage V\ the steady-state tunnelling current den- 
sity at a bias V is I = T/VL 2 S. The result is: 
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where X(9) — (sec 9— l)(sec 2 9 — 2 sec 9— 1) is small for a 
small charge imbalance (9 <§; 1). Thus, the power dissi- 
pation into spinwave pairs is independent of the voltage, 
and so the tunnelling current has a 1/V divergence. This 
should be cut off at low bias when the bias V becomes 
comparable to Vq and perturbation theory breaks down. 
This contribution to the tunnelling current is significant 
at low temperatures: it is our main result, consistent with 
the experimental observation^ of a peak at a bias close 
to the spinwave gap. Its contribution to the tunnelling 
current density is shown in Fig. [21 

We note that the region of the experimental sam- 
ple in which the tunnelling current flows remains an 
open quest ion |4| . Estimating the power dissipated in the 
experiments 2] to be 10~ 16 W, and neglecting renormal- 



isation of the tunnelling (so Asas = Asas), we obtain 
an estimate of 50/zm 2 for this area. This estimate is in- 
creased if we use a renormalised Asas- The result is 
consistent with tunnelling taking place near the contacts 
to the bilayer, rather than over its entire area. In the 
experimental data of 2], the current decays more slowly 
than 1/V away from the resonance, so that the power in- 
creases with increasing applied voltage. We attribute this 
increase to other dissipative channels, due to disorder or 
finite temperature. Processes at higher order in 1/5 also 
affect the current in this way (see {7J) below), but are not 
large enough to explain the differences between theory 
and experiment. 

It should be noted that the dissipation rate © 
is averaged over the period of the Josephson oscilla- 
tions. As with other theories in which the tunnelling is 
treated perturbatively0,||j, we also expect an oscillatory 
a.c. component with frequency eV/Ti. 

At this point, we make contact with previous calcu- 
lations of the tunnelling current in Q . If we ignore the 
weak 0-dependence of the vertex factor, the dynamics de- 
pend only on the azimuthal angle <p of the pseudospin. 
This is the same theory as in jj]. However, that work 
gives a vanishing tunnelling current at zero temperature. 
Our calculation is different from that one since we include 
quantum corrections involving multi-spinwave processes. 

More quantitatively, the current is calculated in 
from Is ~ J d 2 rdt exp(ieVt/Ti — Gq) where the quan- 
tum propagator Gq is, in our notation: 
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We have explicitly included a lattice cutoff, a ~ Ib, that 
reflects the finite bandwidth of the spin waves. We re- 
call that our parameters satisfy (Kv/ peSIb) < 1 to avoid 
quantum disordering. We can therefore perform the in- 
tegral order by order in (fiv / peSo), finding: 

= cM-eVa/hv) DA 2 d 2 ^ 
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where ft — eV/AirSpE and Iq is a modified Bessel func- 
tion. The leading terms in the voltage V are also the 
leading ones in 1/S, justifying the expansion about the 
classical limit. 

In an expansion in eV j PeS, the leading term in eq. Q 
coincides with eq. JSJ if we ignore X(9). In fact, the nth 
term in the expansion arises from the decay of a single 
quantum of the q = mode into n + 2 finite-momentum 
spinwavcs. Thus, the zero-temperature tunnelling cur- 
rent arises from multi-spinwaveprocesses. These contri- 
butions have been ignored in |j| although they are, in 
principle, contained in the framework of that paper. 

We now generalise to the case in which a magnetic field 
is applied in the plane of the bilayer. In that case, the 
tunnelling term in the Hamiltonian acquires a spatial de- 
pendence: Asasw 21 — > AsasI""^ cos(Qx) + m v sin(Qx)], 
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where Q — (eByd/hc) 1 ^ 2 is the wavelength associated 
with the in-plane field B y . The spacing, d, between the 
two electron layers is approximately 30nm. 

The tunnelling vertex now involves a change of mo- 
mentum Q in the x-direction. A new process appears: a 
single q — quantum may decay into a single spinwavc 
(Fig. [2)) . The contribution from this process leads to a S 
function peak in the dissipation. The two-spinwave pro- 
cess contributes to the dissipation rate only for eV > vQ. 
The leading contributions to the current density are 
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where <d(x) is the step function. The peak at eV = vQ 
is asymmetric. The lineshape is controlled by the multi- 
spinwave processes. For eV < vQ, there is no available 
spinwave channel for dissipation and so there is no tun- 
nelling current: at larger voltages, the current decays 
as a power law. This is shown in Fig. |3 in which we 
have broadened the delta function by introducing low- 
momentum scattering with momentum spread gq. This 
may arise from disorder or thermal effects. In the absence 
of these effects, the delta function will be broadened by 
the intrinsic lifetimes of the final states (for example, by 
further tunnelling or spinwave/spinwave scattering). 

The treatment thus far has been valid only for voltages 
greater than Vo = 2A sw /e. We now discuss charging volt- 
ages smaller than this value. In this regime, our thought 
experiment does not result in a steady tunnelling cur- 
rent. Instead, the current and charge imbalance oscillate 
coherently around zero with frequency A sw /7i. We now 
evaluate the rate at which these oscillations decay. 

We treat the parameter 8 perturbatively in this regime. 
In the Holstein-Primakov representation, the quadratic 
part of the Hamiltonian takes the form: 



(9) 



1 /2 

where uo q = [(A SAS + p E -f(q))(D + A SAS + p E j(q))] 
and f3 q annihilates a spinwave with wavevector q. 

The initial state of our thought experiment has a 
charge imbalance. This corresponds to a condensate of 
bosons with q=Q. These bosons have the small concen- 
tration n /S = 9 2 [1 + (£>/A S as)] 1/2 / 47 ^b- The simplest 
dissipative process involves four q—0 quanta annihilating 
to form a pair of spinwaves with finite momenta. We 
defer details to a later paper 01 : the dissipation rate is 
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approximately equal to 3 since As as *C D. This dissipa- 
tion corresponds a decay of the density of q = modes: 
h = —T(9)/ihA sw L 2 . Since the density n is propor- 
tional to 9 2 , we see that d9 2 /dt ~ — 9 8 . The amplitude 
A of the oscillations in the tunnelling current is propor- 
tional to sin#, and so it decays in time as a power law: 
A(t) ~ l/t^e. 

The dissipation is very weak at small initial angles 
partly due to the kinematic constraint for a 4-spinwave 
collision. One might expect stronger dissipation in the 
presence of disorder or inelastic scattering, e.g., with 
phonons or charged quasiparticles. Nevertheless, the 
multi-spinwave processes provide an intrinsic damping 
mechanism that has not previously been identified. 

We note that this treatment does not lead to a state 
with a direct tunnelling current at V < Vq. Therefore 
we cannot address the lineshape of the zero-bias peak in 
dl/dV. It is worth noting that this feature has sharp- 
ened with improved experimental conditions. It may be 
controlled by low-energy states introduced by disorder 
or by effects arising from coupling to external leads. An 
approach similar to [f| may be needed. 

In conclusion, we find that the zero-temperature tun- 
nelling current has a strong 1 /V peak arising from the de- 
coherence of Josephson-like oscillations by the generation 
of electron density fluctuations in the bilayer. We also 
predict a non-trivial lineshape for the dispersing feature 
that appears in the presence of an in-plane field. These 
effects resemble features reported in experiments 0- 

We acknowledge the financial support of EPSRC (in 
part through GR/S61263/01) and of the Royal Society. 



where f(x) = (1- 



l [3-4a;-8a; 2 -.TVl + 16a: + 16x 2 



is 



[9 
[10 

[11 

[12 



S. M. Girvin, A. H. MacDonald and J. P. Eisenstein, in 

Perspectives in Quantum Hall Effects, ed. S. Das Sarma 

and A. Pinczuk (Wiley, New York, 1997). 

I. B. Spielman et at, Phys. Rev. Lett. 84, 5808 (2000); 

Phys. Rev. Lett. 87, 036803 (2001) 

K. Moon et al., Phys. Rev. B 51, 5138 (1995). 

L. Balents and L. Radzihovsky, Phys. Rev. Lett. 86, 1825 

(2001). 

H. A. Fertig and J. P. Straley, Phys. Rev. Lett. 91, 
046806 (2003). 

A. Stern et al, Phys. Rev. Lett. 86, 1829 (2001); M. 
Fogler and F. Wilczek, Phys. Rev. Lett. 86, 1833 (2001); 
Y. N. Joglekar and A. H. MacDonald, Phys. Rev. Lett. 
87, 196802 (2001); 

M. Kellogg et al. Phys. Rev. Lett. 93, 036801 (2004). 
X.-G. Wen and A. Zee, Phys. Rev. B 47, 2265 (1993); 
Z. F. Ezawa and A. Iwazaki, Phys. Rev. B 48, 15189 
(1993). 

H. A. Fertig, Phys. Rev. B 40, 1087 (1989) 
J. P. Eisenstein et al, Phys. Rev. Lett. 68, 1383 (1992); 
S. Q. Murphy et al., Phys. Rev. Lett. 72, 728 (1994). 
A. V. Chubukov, S. Sachdev, and J. Ye, Phys. Rev. B 
49, 11919 (1994). 

R. L. Jack, D. K. K. Lee and N. R. Cooper, in preparation 



